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Causality-based criteria for a negative refractive index must be used with care 
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Using the principle of causality as expressed in the Kramers- Kronig relations, we derive a gener- 
alized criterion for a negative refractive index that admits imperfect transparency at an observation 
frequency oj. It also allows us to relate the global properties of the loss (i.e. its frequency response) 
to its local behaviour at uj. However, causality-based criteria rely the on the group velocity, not the 
Foynting vector. Since the two are not equivalent, we provide some simple examples to compare 
the two criteria. 



This is a slightly longer version of the paper published 
as Phys. Rev. Lett. 101, 167401 (2008); the extra mate- 
rial is largely that removed to fit within the 4-page length 
limit, although fig. 3 was removed as superfluous during 
review. 

Remarkable electromagnetic properties can be seen 
in materials engineered so that the phase velocity of 
electromagnetic-wave propagation opposes the electro- 
magnetic power flow; such materials are often called "left 
handed" (see e.g. 0, i, i, i, i, S, 0, H H ) , but are more 
precisely described as negative phase velocity (NPV) me- 
dia. As might be expected in a rapidly evolving field of 
research, a variety of conditions [lol 11, [l2| for NPV 
have been proposed in the literature. The presence of 
a negative refractive index (NRI) allows for a number 
of intriguing possibilities, e.g.: the creation of a "per- 
fect lens" that produces an undistorted image without 
causing any surface reflections [l3|, the possibility of a 
reversed Casimir force being used to levitate ultrathin 
mirrors [3], the automatic compensation of dissipation 
or dispersion to enhance quantum interference Il5l| . or 
the possibility of "trapped rainbow" light storage [iq . 

The dispersive nature of the effective medium parame- 
ters is exploited in metamaterials to produce an NRl, as 
confirmed through experimental, theoretical, and numer- 
ical studies 0, i, i, i, H, 0, S, [a, [il)- Such metamate- 
rials therefore inherit unavoidable losses on the grounds 
of causality. Since losses can cause a significant drop in 
performance, a key challenge is to successfully compen- 
sate for loss by adding gain; but note that care must be 
taken in theoretical investigations to ensure that the gain 
model is both stable and causa l fisl . ; we also need to 
use the correct NPV criterion [20|. Here we specifically 
address the role of the losses required by causality, by 
considering the famous Kramers-Kronig (KK) relations 



(see e.g. [2l|) which control the relationship between the 



real and imaginary parts of the electric and magnetic ma- 
terial responses (i.e. the permittivity e and permeability 
/i). Such relations can also be established for the square 
of the refractive index — c^e/i, as this quantity in- 
herits the analytical properties of e(cij) and ij.{uj): i.e. it 
lacks singularities in the upper half-plane of complex co 



to place limits on the minimum losses that accom- 
pany NRI for a medium which is perfectly transparent 
at the observation frequency. He concluded that any sig- 
nificant reduction in the losses near the chosen observa- 
tion frequency will also eliminate the NRI. Whether real 
metamaterials can in principle be made with low loss is a 
question of utmost importance in practical metamaterial 
design. Previous work 22| claimed that the answer is 



emphatically negative, but we show here that the answer 
is actually affirmative. 

In this Letter we replace Stockman's zero-loss crite- 
rion with another causality-based criterion, one capable 
of giving useful answers for NPV propagation because it 
admits arbitrary linear optical losses both at and away 
from the observation frequency. Here we assume a ho- 
mogeneous medium with e and /i being effective parame- 
ters obtained for the composite metamaterial by (e.g.) a 
modified S-parameter technique [23, 24]. Such effective 
medium approaches are less reliable in the short wave- 
length (high frequency) regime; but existing analytic at- 
tempts only apply to (at best) thin composite layers [H]. 
The KK relation for n'^{uj) can be written 



ie in" 



Im(n2) 



sds, 



(1) 



where Me() and Im() take the real and imaginary parts; 



thus Im (n^) 



e'n" 



e"/i', where e' and fx' are the real 



parts of e and /x, with the imaginary parts being e" and 
fi"; V takes the Cauchy principal value. Noting that the 
loss in the material is important in the calculation of 
Im {n'^), Stockman transformed Eq. ([T]) into one relating 
the material loss to the presence of NRI, 



= 1 



UpUg 



(s2 _cj2)2 



s^ ds < 0, 



(2) 



by applying the operation L = uj ^d^uj"^ to both sides 
of Eq. 11]). The behavior of several experimental systems 



was claimed in 221 to be consistent with this criterion. 



Here the NRI condition relies on opposed (real valued) 
and n'^{u)) — > 1 as ^ oo (see e.g. 0, IHj). In a re- phase and group velocities, i.e. VpVg < 0. Since the usual 
cent Letter, Stockman ^2^1 adapted the KK relation on acronym NPV is ambiguous, we refer to VpVg < as 
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NPVG (i.e. NPV w.r.t. group velocity); the usual case-^ 
is then NPVE (i.e. NPV w.r.t. energy velocity). 

When VpVg < 0, the integral on the right-hand side 
(RHS) of Eq. ^ must be negative. Stockman there- 
fore concluded that even if the losses vanish at the ob- 
servation frequency, there must still be significant loss 
nearby, otherwise the integral will produce a positive re- 
sult. Consequently, systems with imperfect tuning or 
an insufhciently narrow operating bandwidth would have 
their performance degraded. 

The limitations of Eq. ([2|) are threefold: 

1. Im(n^) and its derivative must be exactly zero at 
the observation frequency uj - otherwise the integral 
diverges, and the constraint becomes uncertain. 

2. It only applies at a particular observation fre- 
quency, despite utilizing the global properties of the 
material response. It can be used to infer the pres- 
ence of nearby loss, but does not indicate whether 
NRI is present there. 

3. The NPVG condition VpVg < is not equivalent to 
the NPVE condition P • fc < 0. 

These limitations make it hard to determine how losses 
might be minimized whilst still maintaining NPV over 
some frequency window. 

We now replace Stockman's criterion with one that 
avoids divergences while allowing for non-zero loss; thus 
removing the first two limitations given above. All nec- 
essary convergence or limiting properties for n'^ can be 
satisfied if e and fi are described by functions of s which 
are both rational and causaP. The third limitation is in- 
trinsic to the approach, but has the advantage that it also 
enables us to evaluate the presence of NRI (or, strictly, 
the presence of NPVG) using causality. Here the group 
velocity Vg amounts to the commonly used d,,,k(uj)\ a,l- 
though imperfect in the case of loss or gain (see e.g. [28|) 
it is that which follows most naturally here. 

Our first step is to integrate the RHS of Eq. ^ by 
parts, but only after splitting it into two pieces covering 
the ranges [0, uj — a) and (cu + a, oo), then taking the limit 
CT — > at the end. After defining Qj = 9^Im (n^(s)), and 
with Qo{s) = Im (ri^(s)) tending to zero fast enough so 



^ The usual NPV condition is P ■ fc < 0, with Poynting vector P 
and wavevector fc. Using the electromagnetic energy density p, 
we can define an energy velocity = P ■ k/ p\k\. This allows us 
to re-express P.fc < as VpV^ < 0. The definition can also be 
extended to embrace moving media, where P is replaced by the 
electromagnetic energy momentum tensor |26ll . 

^ If necessary, inconvenient singularities or divergences in n^(u;) 
can be removed by considering f {uj)'n? {uj) , where / is some ra- 
tional function of u) designed to cancel the pole or to remove the 
divergence |27|| . 



that the s = oo surface term vanishes, we find 

Qo{s) , 
— -s as - 





s2 







ds. (3) 



Since the RHS is independent of Qq, we can now obtain 
a criterion valid where loss is present at the observation 
frequency; and the better behaved integrand means it is 
considerably easier to make inferences about the presence 
of NRI. After applying £, we find that VpVg < requires 



TT< I Qi(s) 

'0 



In 



1 - 



ds. 



(4) 



Before taking a ^ 0, a, second integration by parts 
yields 



TT < 



h{s)ln 

-a 
Q2{S) 



1 - 



ds 



UJ + (7 
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ds 
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(5) 



Again the surface terms at lu — a and uj + a will cancel; 
those at and oo vanish. With z = s/u, we now have 



n < LU 



i{zuj) In |l — z'^l dz 



tanh ^(z) dz. (6) 



The first (Qi) term can again be integrated by parts, 
however the result adds little insight. 

Eq. © is the most general causality-based criterion 
achievable, and is, crucially, not restricted to points of 
perfect transparency. It depends only on how the loss (as 
specified by Qo = Im(n^)) changes with frequency (i.e. 
on its dispersion, as given by Qi and Q2), and not on 
its magnitude. Notably, the sign of Qi (i.e. whether Qo 
is increasing or decreasing with frequency) has a strong 
effect on the presence of NRI; as does the sign of Q2 
(crudely, whether Qo is near a minimum or maximum). 

The non-Qi parts of the integrands (i.e. ln|l — z^| 
and tanh~"'^(z)) are both strongly peaked at s = cj, 
but nevertheless have finite integrals. Using the expan- 
sion Im(n2(s)) ~ Qo{lo) + xQi{uj) + {x^ /2)Q2{uj), for 
X = s/o; — 1 , we can integrate Eq. ^ analytically in an 
attempt to obtain an approximate criterion 

TT < -1.344 ujQi{uj) [I - 2ujQ2{ijj)] - 1.386 uj'^Q2{uj). 

(7) 

Unfortunately this fails to convincingly match Eq. ([6]), 
and the attempt only succeeds in emphasizing that it 
is the global properties of the loss which constrain the 
presence of NRI. Only in Stockman's zero-loss case is a 
simple intuition valid. 
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Using our causality-based criterion in Eq. we can 
now try to infer whether the global properties of Im(n^) 
promote (or hinder) NPVG. Since VpVg and Eq. ^ are 
intimately connected by the KK relations for n?, we used 
this to numerically test the examples below; neverthe- 
less, each expression provides its own unique perspective 
- one local, one global. Since we may not always be able 
to rely on obtaining from a model (as in 23, 23] )j '^^ 
may need to recover it from experimental data. Whilst 
the standard KK relations are prone to generating inac- 
curate reconstructions, approaches such as the Multiply 
Subtractive KK method can resolve this for ma ny p rac- 
tical applications - even nonlinear spectroscopy (29| . 

We now proceed to test our causality-based criterion. 
Since we wish to emphasize general principles, and ensure 
the points we make are clear, we consider simple exam- 
ples with e = jj,, rather than more complicated systems; 
we also normalize with respect to some suitable reference 
frequency. The condition natural to the approach used 



here is the NPVG one (i.e. xg 



Vr,V, 



< 0). This means 



we only need to calculate (and show) one of xg or Eq. ([H]); 
we label the result xg- In contrast, the NPVE condition 
requires that the phase velocity is opposed to the energy 
velocity. This occurs if 



velocity : 

i3,[23| 



XE = e' ImI kl < 0. 



(8) 



The two conditions {xg < and xe < 0) will agree 
if the group velocity Vg and energy velocity ve have the 
same sign. However this only hold in the limit of nearly 
undistorted pulse propagation [30| . i.e. for small disper- 
sion and loss. This is likely to be a poor approximation 
in NRI materials, which by their nature rely on strong di- 
electric or magnetic response. So although our criterion 
in Eq. fH]) can be always used to judge the presence of 
NPVG, and make inferences thereon, this is not strictly 
equivalent to the presence of NPVE. 

Our first exarnple is a simple double-plasmon reso- 
nance, as in e.g. |3l| . setting e and /i according to 
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eo 



Mo 



uj {uj + 17) 



(9) 



A simple test to evaluate the presence of NPVG (and 
at the same time test our generalized causality-based cri- 
terion in Eq. ©), is to increase the losses whilst compar- 
ing it against the NPVE condition. The results can be 
seen in Fig.[l] where x^/"^ is plotted to accommodate the 
vertical range. For sufficiently weak losses (7 ^ uj) the 
criteria agree, with both the xe and xg curves remain- 
ing below zero. However, as the losses get stronger, the 
XG and Xe start to disagree. Nevertheless, we can see 
that in the preferred region of w ~ 1, where e = /i ~ — 1, 
they disagree only for very large losses. Here the xg 
criterion works relatively well because the plasmonic re- 
sponses vary both smoothly and monotonically, hence Vg 
does not change sign and remains in accord with ve- 




FIG. 1: A double (e and /i) plasmon system exhibiting NRI, 
with both plasma frequencies being cup ~ 1.4. It compares 
the (xe) NPVE condition (thick lines), to the (xa) NPVG 
one (thin lines). The results shown are for 7 — 0.02 (solid 
lines), 7 — 0.04 (dashed lines), and 7 — 0.06 (dot-dashed). 



Our next example is again motivated by simplicity, but 
also by the possibility of creating NRI in atomic gases. 
In a gas, it is p ossible to design pumping schemes that 
create gain [3^ , but the freedom to manipulate the 
optical properties relies mainly on the dielectric response 
(e). Here we consider two matched pairs of Lorentz res- 
onances, so that e{Lu) — /i(ti'), and 



eo 



1 



0-1^1 



(72 w| 



«CJ272 



(10) 



We focus on a dominant lossy resonance (cti < 0), 
with a weaker, offset, active resonance ((T2 > 0) providing 
sufhcient gain to induce near-transparency at a chosen 
observation frequency'^, e and /i are chosen equal apart 
from a scale factor eo//io, and are shown on Fig. [2Ia), 
where we see that near transparency has been achieved at 
the cost of increased dispersion, with e' varying strongly 
where uj lu2. Note how the sign of xg swaps back and 
forth according to the gradients of e' and /i', even though 
the values of e' and fj,' themselves change very little: the 
utility of the XG criterion depends entirely on whether 
Vg has the same sign as ve at the frequency of interest. 

The narrowband region of low loss in this system makes 
it ideal for examining our NPVG criterion of Eq. ^ 
in more detail. First, note that there is an asymmetry 
about the loss minimum ~ below, the two contributions to 
Eq. ([6]) reinforce to help satisfy the criterion; above they 
partly cancel, making NRI less likely; this asymmetry is 
visible on Fig. [2lb,c) around uj — 1.05. At the minimum 
itself, we can expect the Qi integral to be small since 
the integrand near uj will be not only small but odd; the 
behaviour will then be dominated by that of Q2 - and 
indeed on Fig. El^c) there is strong qualitative agreement 
between Q2 and xg- The criterion therefore controls the 



^ Note that this has a causal loss profile containing a minimum 
- a situation supposedly excluded by the Eqs. (5,6) and related 
discussion in 1231 • 
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FIG. 2: A system exhibiting narrowband NRI. It combines a 
lossy resonance at tJi = 1 (with 71 = 0.05, ai — —5), and an 
active one at UJ2 — 1.05 (with 72 = 0.01 and (J2 — 1.02). (a) 
The real parts and imaginary parts of e and /i. (b) Com- 
parison of the NPVE (xe) and NPVG (xa) criteria, (c) 
Expanded view around UJ2 = 1.05, showing also the NPVG 
approximation from Eq. (O (labeled xga), and Qi and Q2- 



width of allowed low-loss windows: a narrowband win- 
dow will have a large Q2, so that our criterion will be 
more easily satisfied. This inference is related to Stock- 
man's - it also demands sufficient loss close to the obser- 
vation frequency, but does not require Qo = Qi = 0. 

In conclusion, we have derived a causality-based crite- 
rion for NRI allowing for frequency dependent (disper- 
sive) losses at the observation frequency. Our new crite- 
rion is applicable to any medium with the linear response, 
required by the Kramers-Kronig relations. We inves- 
tigated our causality-based criterion using some simple 
material response models, showing that since the group 
velocity Vg does not always match signs with the energy 
velocity ve, the NPVG and NPVE forms of NRI are not 
equivalent. Since NPVE (i.e. P-k <0) is usually the pre- 
ferred condition for NRI, this difference needs to be taken 
into account before causality-based NRI conditions are 
utilized. Nevertheless, our causality-based NPVG crite- 



rion provides unique insight into how the global response 
of the material affects its local performance. 
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APPENDIX: PREVIOUS COMMENTARY ON 
STOCKMAN RESULT 

Stockman's criterion has been s ubj ect to recent com- 
ment by Mackay and Lakhtakia 



Pendry, Appl. Phys. 
K. Berggren (2007), 
H. Keitel (2007), 



3J|, who constructed 



a dispersion relation designed to provide a counterargu- 



ment. Stockman then disputed [35[ the physical plau- 
sibility of this, and remarked that it did not minimize 
NPV losses very effectively. While the acceptable level 
of losses varies according to the application, the limita- 
tions of Stockman's result make it hard to determine how 
losses might be minimized whilst still maintaining NPV 
over some frequency window. 



APPENDIX: SYMMETRIC LOSS AND GAIN 

Next we combine matched active and passive reso- 
nances (i.e. (Ti = —(72 in Eq. (fTO|) ) offset either side of the 
observation frequency (see Fig. (Ha)), as in e.g. Anant 
et al. [13] ■ This configuration achieves a rather wide 
NPVE bandwidth - if it could be engineered, then it 
would be an attractive system for use. However, just as 
for the offset-gain model above, the NPVG XG criterion 
performs poorly compared to the NPVE one because of 
the way the sign of Vg alternates while following the de- 
tailed behaviour of e' and fi' . 
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FIG. 3: A system exhibiting a wideband and weakly dis- 
persive region of NRI (NPVE) near w ~ 1. It combines a 
lossy resonance at wi = 0.95, ai = 5 and 7 = 0.05, bal- 
anced against an active one at W2 = 1.05, with 02 = —cri and 
72 = 7i- (a) The real parts and imaginary parts of e and /U. 
(b) Comparison of the NPVE and NPVG conditions. 



